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Wilson-like real-space renormalization group and low-energy effective spectrum of 

the XXZ chain in the critical regime 

Kouichi Okunishi 

Department of Physics, Niigata University, Igarashi 2, Niigata 950-2181, Japan. 

We present a novel real-space renormalization group(RG) for the one- dimensional XXZ model 
in the critical regime, reconsidering the role of the cutoff parameter in Wilson's RG for the 
Kondo impurity problem. We then demonstrate the RG calculation for the XXZ chain with 
the free boundary. Comparing the hierarchical structure of the obtained low-energy spectrum 
with the Bethe ansatz result, we find that the proper scaling dimension is reproduced as a fixed 
point of the RG transformation. 
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In the last decade, the density matrix renormaliza- 
tion group(DMRG) has been the most reliable numer- 
ical renormalization group(NRG) for the ground-state of 
the one-dimensional(lD) quantum many body system. 1 
Before DMRG, however, the conventional NRG scheme 
such as a block spin transformation often failed for the 
ID quantum system. The reason for this was analyzed 
carefully by White; 2 the effect of the boundary is signif- 
icant in the ID system, and his DMRG algorithm over- 
come it successfully. On the other hand, Wilson's NRG 
procedure for the Kondo impurity problem 3,4 has been 
highly successful for solving the various impurity prob- 
lems. In Wilson's NRG scheme, the log-discretized energy 
shells for the s-wave electrons are mapped into an effec- 
tive tight binding model with a cutoff parameter, which 
is essentially the ID quantum many body problem with 
a boundary. However, it should be remarked that Wil- 
son's NRG is efficient for the gapless system, in contrast 
to the DMRG which prefers the gapful system. 

In this paper, we focus on the NRG for the gapless ID 
quantum spin system; Wilson's NRG for the impurity 
problem is reconsidered in the context of the ID quantum 
spin chain. The essence of the idea is very simple; remove 
the impurity site and add the interacting spins instead of 
the free electrons in the impurity problem. Such an ap- 
proach without the cutoff parameter was tested by Xiang 
et al in the almost the same timing as the DMRG, but 
it seems more suitable for the gapful spin chain. 5 Also, 
the infinite-system-size DMRG still cannot overcome the 
critical fluctuation in the ID critical system in the bulk 
limit; An essential point in the present approach for the 
critical system is that we introduce the cutoff parameter 
controlling the energy scale of the system, which plays 
an important role in the RG transformation. 

In the following, we briefly formulate the Wilson like 
NRG for the XXZ spin chain, where we emphasize the 
role of the cutoff parameter. Numerical calculations are 
actually performed for the XXZ chain in the critical re- 
gion, and the obtained results are analyzed on the basis 
of the Bethe ansatz/CFT solution. We then find that the 
low-energy excitation spectrum of the XXZ spin chain is 
successfully reproduced as a fixed point of the RG trans- 
formation. 



In this paper, we deal with the 5 = 1/2 XXZ chain in 
the critical regime for simplicity, but the formulation for 
the general case is straightforward. We write the local 
Hamiltonian of the XXZ chain as 



U QX QX 

ll>n,n+l — d n °n+l 



CV CV 
a n J n+l 



(1) 



where S is the 5 = 1/2 spin matrices and < A < 1 is 
assumed. We then consider the Hamiltonian of N spins 
with the cutoff A as 



H N (A) 



N-l 

53 A N - n - i h 

n=l 



n,n+l ; 



(2) 



for which the free-boundary condition is basically as- 
sumed. If A = 1, eq. (2) becomes the uniform XXZ chain 
with the free-boundary condition. In Wilson's original 
NRG for the Kondo problem, there is an impurity at the 
n = site, and /i n ,n+i should be the tight binding free 
electrons. In the present case, n — is the empty site 
and the interacting spins are adopted as h n , n +i. 

In order to treat matrices having different sizes in the 
RG transformation, we use the notation for a m x m 
matrix X: 



X* = I ® X, 



(3) 



where / is the 2x2 identity matrix for the space of the 
spin added, and thus, X* becomes a 2m x 2m matrix. 
The recursion relation between the N- and N + 1-spin 
systems is given by 



H N+1 (A) = AH* N (A) + h NiN 



(4) 



In this recursion relation, the smallest energy scale is 
fixed to be in the order of unity. If A = 1, we do not 
touch the energy scale of the system, leading the simple 
recursion relation for extending the system size. For A > 
1, a new spin having a smaller energy scale is added to 
the bulk part of the system. Thus the outer spins(the 
smaller index of n) has the higher energy scale. 6 

Now we convert the basis so as to diagonalize Hjy(A), 



H N (A)U N = U n uj n , 



(5) 



where u) are the eigenvalues of the Hamiltonian and U 
are the corresponding eigenvectors. By maintaining the 
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lower energy states on the basis of (5), we obtain the 
recursion relation for the renormalized Hamiltonian: 



H N+1 (A) = Al 



'N 



Am 



(6) 



where i?jv+i(A) = llj^ H N+1 (A)U'^. By diagonalizing 
the Hamiltonian of N sites, we can then increase the sys- 
tem size through eq. (6) with maintaining lower energy 
states. After a sufficient number of iterations starting 
with the initial Hamiltonian of N = 2, the matrices con- 
verge to those in the N — > oo limit. However, we should 
note that the energy scale of the system is reduced by 
the cutoff parameter A, which is a crucial point in the 
RG for the critical system. 

-XY model (A = 0)- Let us demonstrate the NRG 
for the XXZ chain in the XY-like anisotropic region(0 < 
A < 1), for which the ground-state has the gapless ex- 
citations. In particular, we first analyze the free fermion 
case(A = 0). Figure 1 shows TV-dependence of the low- 
est eight eigenvalues AE\ of the XY model with no 
cutoff(A = 1) and with cutoff A = 1.05. In the figures, 
AE\ is measured from the "ground-state" energy. We 
can see the oscillation with respect to N =even or odd, 
for which the total S z of the system takes an integer 
or half integer alternately. Although there are the eight 
eigenvalues in the figure, we can see only 5 lines because 
of the degeneracy. 

As is shown in Fig. 1(a), AEa for A = 1 shows clear 
1/N dependence up to JV ~ 50, which corresponds to 
the size scaling of the finite-size energy gap. In principle, 
this implies that the low-energy eigenvalues accumulate 
to AE = in the N — ► oo limit. However, the accuracy 
for N > 50 is not maintained because of the truncation 
of the bases, and thus the lines in Fig. 1(a) show the 
winding behavior. We should note that the range of the 
1/N dependence is limited to up to N ~ 50 even with 
m = 3000. 

In Fig. 1(b), we show the result obtained with the fi- 
nite cutoff A = 1.05, which converges well with respect to 
m(~ 3000). In the figure, we can see that the eigenvalues 
rapidly approach constant values. The regularization due 
to A leads to the NRG iteration converging to the fixed 
point. Since the energy scale of the system is reduced by 
A in each iteration step, this fixed point of the spectrum 
implies that the system is scale invariant. More inter- 
estingly, the energy levels at the fixed point show very 
regular behavior; we can see that the level intervals are 
constant. Here, it should be recalled that the smallest 
energy scale at N = 200 is of order of A -200 ~6x 10 -5 
compared with the largest exchange coupling. 

Figure 2 shows the low-energy eigenvalue distribution 
after sufficient iterations^ = 200). We can see a clear 
stairway-like structure having the almost the same level 
space, which suggests that the fixed point of the free 
fermion model is properly achieved. Here we note that 
this stairway structure is confirmed to be basically the 
same as that for A = 1.03. In order to analyze the 
structure, we invoke the finite size spectrum of the XY 
model, which is obtained by the Jordan- Wigner transfor- 
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Fig. 1. Size dependence of the lowest 8 eigenvalues of the Hamil- 
tonian for the XY model: (a) A = 1 and (b) A = 1.05. The 
ground-state corresponds to the horizontal axis. The inset is a 
log-log plot of the main panel of (a). 



mation, 



Hl — E 



7TV 

T 



n=0 n=Q 

'(7) 

where v is the spin wave velocity, c n and d n are the 
fermion operators, and Eq is the ground-state energy in- 
cluding the 1/L correction. Note that the total S z of the 
system is represented as M — J2 n =o c n c « — d^dn. 

The energy scale in the finite-size spectrum is speci- 
fied by the inverse system size 1/L, while the energy scale 
at the fixed point of the NRG is governed by the cutoff 
log A. However, the degeneracy of the low-energy excita- 
tion in the scaling regime should refer only the scaling 
dimension. Applying Hl to the low-energy states labeled 
by certain fermion numbers, we can identify the degen- 
erating structure of eq. (7) (See also eq. (9) with £ = 1) 



degeneracy # = {1, 2, 1, 2,4,4, 5, 6, • • •}, 



(8) 



which is successfully reproduced in Fig. 2. In addition, 
it can be easily verified that the contribution from the 
sectors labeled by the magnetization quantum number 
M is also consistent with the exact spectrum. 

-XXZ model- We further analyze the XXZ model, 
which contains the nontrivial result of the interaction 
effect. The exact finite-size spectrum of the XXZ chain 
with the open-boundary condition in the critical regime 
was calculated in refs. 8 ~ 10 In particular, the excitation 
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Fig. 2. Low-energy spectrum of the XY model at the fixed point 
for A = 1.05. 



spectrum for the free-boundary case is given by 

AEl(M,N) = —Q(M,N), (9) 

where v is the spin wave velocity and the scaling dimen- 
sion is defined as 



Q(M,N)=^+N. 



(10) 



The quantum number M is an integer corresponding to 
the total S z of the system, and N is a non-negative in- 
teger. In addition, the degeneracy for each N is given by 
Euler's partition number r](N) which is defined by 



n=l 



n=0 



(11) 



In eq. (10), the dressed charge £ can be exactly calculated 
by the Bethe ansatz for the periodic system, 11 

-1/2 



2(1 



arccos A 



(12) 



which represents the nontrivial consequence of the in- 
teraction. Note that the dressed charge of the XY 
model(A = 0) is given by £ = 1, and then eq.(9) can 
be reduced to eq. (7). 
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Fig. 3. The low-energy spectrum of the XXZ model of A = 0.5 
at the fixed point of A = 1.03. 



Figure 3 shows the lowest 50 eigenvalues of the XXZ 
model with A = 1/2, where the magnetization M of the 
system is indicated by the distinct symbols. This spec- 
trum is obtained with the cutoff parameter A = 1.03, af- 
ter 300 NRG iterations with the number of retained bases 
m = 3000. Although the degeneracy of the eigenvalues is 
weakly disturbed by the effect of the finite cutoff, we can 
clearly recognize the stairway structure. In Fig. 3, an 
important difference from the free-fermion case is that 
the level space is not uniform, which is attributed to the 
dressed charge. 

For the case of A = 1/2, the exact dressed charge be- 
comes £ = V3/2. Substituting £ = V3/2 into eq. (9), 
we can resolve the stairway structure of the spectrum 
obtained by NRG; The tower structure of the weight 
Q(M,N) for A = 1/2 is illustrated in Fig. 4. The num- 
bers in the circles indicate the degeneracy correspond- 
ing to r)(N). Moreover, the lowest energy levels in the 
M = ±1 and ±2 sectors are located at Q = 2/3 and 8/3 
respectively. Indeed, the numerical data in Fig. 3 yields 
the ratio AE A (0, 1)/A.E A (1, Q) = 0.68, which is in good 
agreement with the exact value Q = 2/3. By counting 
numbers of the degeneracy for each Q level, we obtain 



degeneracy # = {1,2, 1,2,2,6,3,8,5, • • •}. 



(13) 



These degeneracy numbers are clearly consistent with 
the stairway structure in Fig. 3. Moreover, in Fig. 3, the 
eigenvalues in the M = sector(open circle) clearly show 
(quasi)degeneracy corresponding to rj(N). The M = ±1 
and ±2 sectors also exhibit the same hierarchical struc- 
ture, starting from Q = 2/3 and 8/3 respectively. 
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Fig. 4. Tower structure of the scaling dimension Q(M, N) for 
A = 1/2. The numbers in the circle represents the degeneracy 
numbers r](N) for a given set of quantum numbers (M, N). 



We now turn to an application of the NRG. The ratio 
of the first to the second excitation energy corresponds 
to the dressed charge: 



'!A£a(0,1) 
2 A£ A (1,0)' 



(14) 
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Fig. 5. Dressed charge £ estimated from the NRG spectrum. The 
solid line indicates the exact solution of eq. (12). 



This relation enables us to estimate the dressed charge 
from the NRG spectrum. The estimated £ by NRG up to 
m = 4000 and A = 1.02 is summarized in Fig. 5, where 
the exact solution is also indicated by the solid line. The 
NRG result is in good agreement with the exact solu- 
tion for A < 0.7, where we have clearly confirmed the 
stairway structure of the spectrum. For 0.7 < A < 1, 
however, the open circles slightly deviate from the exact 
line. In this region, we can see that the stairway structure 
in the NRG spectrum is disturbed significantly within 
m = 4000. For a more accurate estimation of £, a A 
closer to unity and a "balanced" increase in m are re- 
quested. However, it should be noted that, at A = 1, 
the degeneracy due to the SU(2) multiplcts properly 
emerges in the spectrum, and thus the correct value of 
£ = l/\/2(— 0.7071) is obtained. This is because the 
SU(2) symmetry is well maintained in the NRG itera- 
tions, although the cutoff A is introduced. 

To summarize, we have formulated the Wilson-like 
real-space RG for the ID XXZ model by introducing the 
cutoff parameter A. We have demonstrated the NRG cal- 
culation for the XXZ chain for < A < 1. Then it was 
confirmed that the hierarchical degenerating structure 
of the low-energy spectrum is consistent with the exact 
scaling dimension which is calculated by the Bethe ansatz 
for the free boundary. Moreover, the estimated dressed 
charge is in good agreement with the exact value. These 
imply that Tomonaga-Luttinger liquid universality can 
be reproduced as the fixed point of the real-space RG 
transformation. In Wilson's original NRG for the impu- 
rity problem, the cutoff A comes from the discretization 
of the Fermi sea in the vicinity of the Fermi surface. On 
the other hand, the cutoff in the present NRG for the 
XXZ model was directly introduced by hand. In both 
cases, the cutoff parameter successfully controls the en- 
ergy scale of the system. Nevertheless, we should note 
that the cutoff also induces an undesirable perturbation 
in the degenerating structure, in addition to the desired 
regularization of the gaplcss spectrum. In order to ex- 
tract the correct scaling dimension, this perturbation 
must be irrelevant compared with the level spaces of the 
desired hierarchical structure. In principle, we can take a 
sufficiently small A(— > 1), but the systematical A depen- 
dence is still unknown. In this sense, further investigation 



of the A dependence is clearly needed. 

As can be expected, it is also straightforward to apply 
the present NRG to the gapful spin chain, and we have 
actually examined the XXZ chain with A > 1. However, 
the obtained spectrum does not have a regular distribu- 
tion, and thus we can not resolve any algebraic structure 
at the present stage. This may be because the regular- 
ization by A competes with the intrinsic energy gap in 
the original XXZ chain. Whether we can extract some 
information on the gapful system using the present NRG 
is an interesting future problem. 

Finally, we make a comment on the infinite-system- 
size method of the DMRG. Of course, the most signifi- 
cant difference of the DMRG from Wilson's NRG is that 
the RG transformation is constructed through the diag- 
onalization of the reduced density matrix. However, an- 
other important aspect particularly for the critical sys- 
tem is that the DMRG does not seem to include any scale 
transformation in its algorithm(the system size increases 
linearly in DMRG). In other words, there is no explicit 
cutoff parameter that controls the energy scale of the 
system. In this sense, we may say that the fixed point of 
the infinite-system-size DMRG is the self-consistent solu- 
tion of a certain matrix product type wavefunction, apart 
from the usefulness of the DMRG. In order to resolve the 
critical behavior in the DMRG, it may be a hopeful way 
that the direct RG analysis of the reduced density matrix 
or the corner transfer matrix, 12 which show the similar 
degenerating structure in their eigenvalue spectrum. 12-15 
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